Experimental Determination of Entanglement for Arbitrary Pure States 
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We present a way of experimentally determining the concurrence in terms of the expectation 
values of local observables for arbitrary multipartite pure states. In stead of the joint measurements 
on two copies of a state in the experiment for two-qubit systems [S. P. Walborn et al. Nature 
(London)440, 20(2006)], we only need one copy of the state in every measurement for any arbitrary 
dimensional multipartite systems, avoiding the preparation of twin states or the imperfect copy of 
the state. 
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I. INTRODUCTION 

Quantum entanglement is one of the most fascinat- 
ing features of quantum theory [![. To characterize 
and quantify the entanglement some well defined mea- 
sures such as entanglement of formation (EOF) 0, Q 
and concurrence [3, [5[ have been used. For given two- 
quibt or some special symmetric states explicit ana- 
lytic formulas for EOF and concurrence have been found 
m S EL S E3, EH ■ For arbitrary given states the en- 
tanglement can be estimated by analytic lower bounds 

Gail M EE M E3 EE M M M . 

Nevertheless, for unknown quantum states, to charac- 
terize the entanglement one needs experimental measure- 
ments. An important approach to detect entanglement 
is the Bell-type inequalities [H, [H H, HE SH . For in- 
stance, Gisin proved that all two-qubit pure entangled 
states violate the Clauser-Horne-Shimony-Holt (CHSH) 
inequality |24( and Chen et al. presented a Bell-type in- 
equality that would be violated by all three-qubit pure 
entangled states (25[. For general mixed two-qubit states, 
Yu et al. [26[ proposed a Bell-type inequality that gives 
a sufficient and necessary criterion for separability. An- 
other experimentally plausible approach is the entangle- 
ment witness (27|, which could also be used to detect 
certain kinds of entangled states with the present tech- 
nology m. 

However, to detect the entanglement in terms of Bell- 
type inequalities one needs expectation values of two or 
more observables (two or more setting measurements) per 
party. And one has to do infinitely many dichotomic 
measurements theoretically. Moreover untill now we still 
have no necessary and sufficient Bell inequalities to detect 
the entanglement for general multiqubit systems. Certain 
entanglement witness works only for some special states. 

In fact the concurrence is defined for both bipartite 
and multipartite states and gives rise to not only the 
separability, but also the degree of entanglement (at least 
for arbitrary dimensional bipartite states). The problem 
is how to use this measure to determine the entanglement 
for unknown quantum states experimentally. 

In [29[ Mintert et al. proposed a method to measure 



the concurrence directly by using joint measurements on 
two copies of a pure state. Then Walborn et al. presented 
an experimental determination of concurrence for two- 
qubit states [30l . [3lj , where only one-setting measurement 
is needed, but two copies of the state have to be prepared 
in every measurement. 

In this paper, we give a way of experimental determina- 
tion of concurrence for two-qubit and multi-qubit states, 
such that only one-copy of the state is needed in every 
measurement. To determine the concurrence of the two- 
qubit state used in (30l . [3lj , also one-setting measurement 
is needed, which avoids the preparation of the twin states 
or the imperfect copy of the unknown state, and the ex- 
perimental difficulty is dramatically reduced. As exam- 
ples general two-qubit and three-qubit systems, and gen- 
eralized multiqubit GHZ states = ao|0 • • • 0) + |1 • • • 1) 
are investigated explicitly. The results are then gener- 
alized to the case of arbitrary dimensional multipartite 
pure states. 



II. CONCURRENCE FOR N-QUBIT SYSTEM 



For a TV-partite M dimensional pure state 



■■■ , i N =0 ••• . ' j W/j «u, 

concurrence is given by $2] , 



€ C, the 



(1) 



where the summation goes over all 2 N — 2 subsets of the 
N subsystems, pi is the corresponding reduced density 
matrix with respect to the ith bipartite decomposition. 
Up to a constant factor, C(\tp)) can also be written as 



cm) = 



A I 



. 2J J! \ a af3a a >/3' - a af 3'a a >p\ 2 , (2) 

\ P {a, a', 13, 13'} 



where a and a' (j3 and j3') are subsets of the subindices of 
a, associated with the same sub-Hilbert spaces but with 



2 



different summing indices, a (a 1 ) and (3 (/?') span the 
whole space of a given subindex of a. ^2 p stands for the 
summation over all possible combinations of the indices 
of a and (3. 

Our main aim is to re-express the concurrence in 
terms of the expectation values of local observables. 
We first give a general proof that this can be always 
done: the squared concurrence of 7V-qubit pure state \tp), 
C 2 (| - 0)), can be expressed by the real linear summation 
of M \ib)(ib\o4, a,-, • • -en 
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where the coefficients x 
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(3) 

j N are real, ctq is 



the 2x2 identity matrix, a\ 
1 
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. 0"2 = ( 
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and cr 3 = ( q ^ ) are the Pauli matrices. 



We only need to show that each term in the squared 
C(\tp)) of can be written in the form of the right hand 
side of ©. Note that 



\Q>apQ>a'P' ~ O-aP'O-a'p] 2 = |aa/9Oa'0' I + \ a aj3' O-a' pY 

- a* a p a *a'i3' a ai3'aa'fi - a al3 a a < p, a* p , a* J(3 . 



(4) 



Set A( n ) = \a/3}(al3\, = \a'/3')(a'/3'\, A( 12 ) = 

\a(3'){a(3'\, A^ = |a'/3)(a'/?|, then 



Kpaatpl* + \a af3 ,a a ,tf = |^}(^ (21) |V> 

-h<-0l-4 (12) |-0> <-0|v4f 22 > |^>- 

(5) 

AW , i,j = 1, 2, obviously has the form AW') = | <g> 

■• • ® | * at) (*iv | , where ii, • • • , ja? take value or 1. As 
= 2(00 + 0-3) and |1)(1| = i(cr - o- 3 ), we have 
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for some real coefficients Xi lt ... ! i N ,j 1 ,--- ,j N . -^(\a(3)(a(3'\ — \a(3')(af3\), A^ = A^(\a' (3'\ — 

Denote further A^ = -l={\ a {3)(a{3'\ + \ a p')(a0\), \a' (3 1 ) {a! (3\) , then 
A^ = i(\a'/3){a>l3'\ + \o/p){c/p\), A^ = 



- a* ap a* a , p ,a ap ,a a ,p - a a pa a , p,a* ap ,a* a , p = -((^ 13 ) \tf,) ^\A^) |^) + ^\ A ^ \t/,) {^\A^ |$). (6) 



It is clear that |a)(a|, \a')(a'\, and \0'){/3\ are 

tensor products of|0)(0|, |1)(1|, |0)(1| = -ia 2 ) and 
|1)(0| = 3(0-1 + icr 2 ). Without loss of generality we as- 
sume A( 13 ) = \ii){h\ ® ••• <8> |i a )(i a | <8> (|i 8 +i)0Wi| ® 
• • • <8> |*jv) 0"jv| + |is+i)<i s +i| ® ••• ® |jjv)(«jv|), where 
1 < s < N, ik,jk take values or 1 and it- ^ jk for 
each s + 1 < k < N. The part <8> ••• <8 j is 

the real linear summation of tensor products of co and 
cr 3 . While the rest part T = |z s+ i)(j s+ i|(g)- • -(g) |*jv)C?Jvl + 



I 

\is+i){i s +i \ ® ••• ® |i'iv)(iiv| can be written as 

JV JV 

23^7 (^+^(- i r^)+ 2 ^ I (8) (crr+^-l) 1 ^^), 

Z=s+1 l=s+l 

where pi takes values or 1 for each I. T is further of the 
form 



N-s N- 

2^ 
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< l m < N — s — I. If N — s — I is even, then i N s 1 is real and each factor of a aj is real. If N — s — Z 
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is odd, then (-1)'- + (-\)N-s-i-i m = Q Hence 4(13) ig 
the real linear summation of the tensor products of <7j, 
< i < 3. 

Similarly one can show that A^ 14 \ A^ 23 ^ and A^ 24 ^ are 
real linear summation of tensor products of <7i, < i < 3. 
Thus -a* al3 a* alf3l a a f}:a al i3-a a i3aa'i3'a* al3l a* a ,p and Eq. (j4|) 
can be expressed in the form of real linear summation of 

(V>K • • • ct ijv I V>) (V'kji °"j2 • ■ • IV>) ■ 

Therefore the squared concurrence of iV-qubit pure 
states C 2 (\ip)) can be expressed as the expectation values 
of tensor products of er^ (0 < i < 3), though such expres- 



sions may be not unique (From ([5]) and ([6]) one sees that 
it is possible to find an expression that is invariant under 
the permutations of the N observables) . 

A. Concurrence for two-qubit system 

For any two-qubit state \ip) = aoo|00) + a i|01) + 
aio|10) + an|ll), |a o| 2 + Ki| 2 + |«io| 2 + |aoi| 2 - 1, 

C 2 = 4|a 00 aii - a iaio| 2 , (7) 

which can be expressed as 



C 2 = -(1 + ((T3<t 3 ) 2 - (a 3 <J ) 2 - (<to<t 3 ) 2 - (tT (7i) 2 + (0-3CT1) 2 - (<7 Q a 2 ) 2 + (cr 3 CT2> 2 )- (8) 



Therefore for experimental determination of the concur- 
rence, one only needs to measure (a 3 a 3 ), {<J 3 a\) and 
(0302) respectively. One may also find alternative ex- 
pressions with symmetry under the exchange of the two 
qubits [H]. 

For states in Schmidt decomposition, = ao|00) + 
ai|ll), I a-o 1 2 + |ai| 2 = 1, we have 



C 2 



;(1 + ((T 3 0'3> 2 - (0W3> 2 - (o-30-o) 2 )- 



(9) 



In this case experimentally we only need to measure 

(0303), or simply count the probability P(++), P{ ) 

of projections | ++)(+ + |, | >( | with |+> = 

^(|0) + |1)) and I-) = ^(|0) - |1)) respectively, since 

C 2 = 16P(++)P( — ). For the state a|01) + (3\10) used 
in , it is also true that only one-setting measurement 
is needed. But here we only need one copy of the state in 
every measurement, while in [30] joint measurements on 
two copies of the state are needed in every measurement. 

For small deviation \tp'} = y/1 — e\ip) + \ft\4>) from an 
ideal pure state \ip) due to imperfect preparation, where 
e € 1R and \4>) is an arbitrary pure state, our protocol 
shows that the concurrence obtained from the experi- 
ment is exactly the one of \tp'}- Hence if the parameter 
e is small enough, the difference of the concurrence be- 
tween and would be small enough. For a two- 



crystal type-I down-conversion source, with improper 
spatial mode matching and spectral filtering, the imper- 
fect preparation procedure could result in mixed states, 
p = (1 - e)|V)(^| + e(\a\ 2 \HH)(HH\ + \(3\ 2 \VV)(VV\) 
instead of the ideal pure state — a\HH) + fi\VV), 
where H and V stand for horizontal and vertical linear 
polarizations respectively. That is, the phase coherence 
between \HH) and \VV) is reduced by 1 — e. Therefore 
the actual concurrence of p is smaller than that of \ip), 
C(p) = (1 — e)\a/3\ [31I ? ]. If we still measure the state 
according to ([8]) or (JOj) , we have C(p) = |a/3|. Thus the 
relative error due to mixing is linear in e. 

Remark In principle one can always use tomography 
to reconstruct the unknown state. However it requires a 
large number of measurements. In particular one needs 
S^-setting measurements to reconstruct an arbitrary N- 
qubit density matrix. To obtain all 16 expectation val- 
ues of the two-qubit density matrix, nine-setting mea- 
surements have to be used 33j. From ([8]) we only need 



three-setting measurements to quantify the entanglement 
of the state, which is much simpler than tomography ap- 
proach. 

B. Concurrence for three-qubit system 

For any pure three-qubit state \ifj) = Yli j k=o a ijk\ijk), 
the squared concurrence is of the form, 



C 2 = 4(|aoooam — &ooiaiio| 2 + |aoooani — aoioaioi| 2 + |aoooani ~ aoiiaioo| 2 + |aooiano ~ aoioaioi| 2 (10) 
+ |aooi a no — aon a ioo| 2 + |aoio a ioi ~ aon a ioo| 2 ) + 8(|a ooaoii — a 00 iaoio| 2 + |aooo a ioi — aooi a ioo| 2 
+ l a ooo a no — aoio^iool 2 + l^ooiain — aon a ioi| 2 + l a oio a m ~ aon a no| 2 + l a ioo a m — a ioi a no| 2 )- 

I 



Up to a constant factor, C 2 can be expressed as 
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C 2 = -(9 - 5(ct ct 3 ct ) 2 - 5(ct ^o^3) 2 - 5(a 3 CT ffo> 2 + + (^3^3^o> 2 + + 3(ct 3 o-3CT3> 2 (H) 

-3(ct cto(Ti) 2 - 3(er ericr ) 2 - 3(criCT cro) 2 - (ctoo- 3 cti) 2 - (cio- o"3) 2 - (ctso'ictq) 2 + 3(a aia 3 ) 2 
+3(ct 3 ct (Ti) 2 + 3(erier 3 er ) 2 + (c 3 cr 3 cri) 2 + (ct 3 ctict 3 ) 2 + (cri<T 3 er 3 ) 2 - 3(a aQa 2 ) 2 - 3(cr cr 2 cr ) 2 
-3(a- 2 croC7o) - (<Toa 3 a 2 ) - (c 2 crocr 3 ) - (CT3CT2C0) 4- 3 (0-002 03) 
+3(a 3 a (T2) 2 + 3(o- 2 c 3 <to) 2 + (03C3C2) 2 + (cr 3 cr 2 cr 3 ) 2 + (a 2 a 3 a 3 ) 2 ). 



which is invariant under the permutations of the three 
qubits. For experimental determination of the concur- 
rence for arbitrary three-qubit states, seven quantities are 
needed to be measured: (cr 3 a 3 cr 3 ), (<7 3 cr 3 cri), (cr 3 aia 3 ), 
(aia 3 a 3 ), (cr 3 cr 3 (7 2 ), (cr 3 c 2 cr 3 ), (cr 2 cr 3 cr 3 ). 

In particular for the three-qubit generalized GHZ state, 
\ip) = a |000) + aijlll), |a | 2 + |ai | 2 = 1, and the 
generalized W state \ip) = a |001) + a x |010) + o 2 |100), 
l a o| 2 + | fl i| 2 + l a 2| 2 = 1, their squared concurrence are 
12|a ai| 2 and 8(|aoai| 2 + |a a 2 | 2 + |aia 2 | 2 ) respectively. 
The concurrence of both generalized GHZ states and gen- 
eralized W states can be measured according to the fol- 
lowing formula: 



C 2 = i(9-5(a a 3( 7 ) 2 
+ ((7Q(7 3 a 3 ) 2 + (a 3 a 3 a ) 



5(a a a 3 ) 2 - 5(a 3 a a ) 2 
4- (a 3 a cr 3 } 2 + 3(a 3 a 3 a 3 } 2 ). 

(12) 



(112)) shows that for experimental determination of en- 
tanglement for these states, one needs only one-setting 
measurement, (a 3 a 3 a 3 ). 



Similar results can be obtained for multiqubit systems 
such as iV-qubit generalized GHZ state \ip) = ao|0 • • • 0) 4- 
ai 1 1 • • • 1) , |ao| 2 + I &i 1 2 = 1, or iV-qubit generalized W 



state \ip) = Oo|0 

|2 1 |„. 12 , . 



•01) + ai|0--- 10) + - • - + 0^-1110 ■••0), 
, a or 4- I ail" + ■ • ■ 4- |aAr-i| 2 = 1. For instance for the 
generalized GHZ state, the concurrence is |aoai| up to 
a constant. Its squared concurrence can be expressed as 
follows: 



Kfe,fe'<JV 



C z = 1 



E 



r (»l— ik)\/ rr Ui 

'3 I \ a 3 



k.k' 



is even 



■3k' 



1<1,1'<N 

E < 



■3l')\ 



(13) 



w 



odd 



where (a 3 n ) denotes the expectation value of the lo- 
cal operators such that the ii-th, • • • , ik-th are a 3 oper- 
ators and the rest are identities. 



III. CONCURRENCE FOR TV-PARTITE 
M-DIMENSIONAL SYSTEM 

Besides qubit systems, our approach can be also used 
for arbitrary M-dimensional cases. In stead of the Pauli 
operators, one can use the SU(M) generators as observ- 
ablcs: 

M-i 



s-1 

X s =^2\j)(j\-s\s)(s\, 1<8<M-1, 

1 

A 8 = |iX*l + l*>0'l> s = M,-,-(M + 2)(M-l), 

X s = -i(\j)(k\-\k)(j\), «=i(M + l)M,. - ,M 2 -l, 
where 0<j<fc<M-l. Note that 
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1 



Ao 



1 



M M(M - 1) 



Am-i + 



1 



Am-i 



(M - l)(M-2) 
1 

(M — 1)(M - 2) 



Am-2 + ■ 
Am-2 + ' 



3^2 2 Al ' 



|M - 2)(M - 21 = — A + — r 

\M-l){M-l\ = ^\ -±\ M -i, 



1 



M-l 



M - 1 



VAf-2) 



and for < j < A; < M - 1, = |(A S + iA,/), 

|fc)(j| = I(A s -iA s /) for someM < s < ±(M + 2)(M-1) 
and |(M + 1)M < s' < A/ 2 - 1. Similar to the proof 
of iV-qubit system, it is direct to show that the squared 



concurrence of the A^-partite M-dimensional pure state 
\tp) can be expressed in terms of real linear summation 
of (^\\ ll X t2 ■••AjJ'0)(V>|A.,- 1 A i2 ■ • • \jx\ip): 



M -1 

ii, ••• , in, ji, ■■■ . jn=0 



lAiiAij •••A ijv |i/')(V'|Aj 1 A i2 •••A 



(14) 



where a;^, ... . i N> j u ... t j N are real. 

IV. CONCLUSIONS 

We have proposed a method for experimentally de- 
termining the concurrence in terms of the expectation 
value of local observables, which gives not only sufficient 
and necessary conditions for separability of the quan- 
tum states, but also the relative degree of entanglement. 
Moreover unlike the case of Bell-type inequalities where 
measurements are needed with respect to infinitely many 



observables, we need only mean value of a few observ- 
ables. And in stead of the joint measurement on two 
copies of the state needed in the experiment (2f| [53, HH 
for two-qubit states, we need only the usual measure- 
ments on one copy of the state in every measurement 
for any arbitrary dimensional multipartite states, which 
dramatically simplifies the experiment and reduces the 
error rates and the imperfectness in the preparation of 
the states. Compared with entanglement witnesses, for 
which some a priori knowledge about the states under 
investigation is needed, we do not need any information 
before measuring the state in experiment. 
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